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Abstract. In view of the usefulness and importance of kinetic equations in cer- 
tain physical problems, the authors derive the explicit solution of a fractional 
kinetic equation of general character, that unifies and extends earlier results. 
Further, an alternative shorter method based on a result developed by the au- 
thors is given to derive the solution of a fractional diffusion equation. 

1 Introduction 

Fractional reaction/diffusion equations involve fractional derivatives with re- 
spect to time and space and are studied to describe anomalous reaction/diffusion 
of dynamic systems with chaotic motion. Fractional kinetic equation for Hamil- 
tonian chaos is discussed by Zaslavsky (1994). Solutions and applications of cer- 
tain kinetic equations are studied by Saichev and Zaslavsky (1997). Solutions of 
a fractional kinetic equation is investigated by Haubold and Mathai (2000) for 
a simple production-destruction mechanism. This equation was generalized by 
Saxena, Mathai, and Haubold (2002). In recent articles, Saxena, Mathai, and 
Haubold (2002, 2004a, 2004b) discussed the solution of a number of generalized 
fractional kinetic equations . In the present article, we investigate the solution 
of a unified fractional kinetic equation, which provides unification and extension 
of results on fractional kinetic equations given earlier by Haubold and Mathai 
(2000) and Saxena, Mathai, and Haubold (2002, 2004a). We also present the 
solution of a fractional integral equation discussed by Miller and Ross (1993). 
Further, an alternative proof of the solution of a fractional diffusion equation 



1 



given earlier by Kochubei (1990) is investigated, which is based upon a result 
given by Saxena, Mathai, and Haubold (2006). Most of the results are obtained 
in terms of generalized Mittag-LefHer functions in elegant and compact forms, 
which are suitable for numerical computation. 

The paper is organized as follows. Section 2 contains the solution of a unified 
fractional kinetic equation while Section 3 considers special cases of the equation. 
A shorter alternative method for the solution of a diffusion equation discussed 
earlier by Kochubei (1990) is presented in Section 4. A series representation 
and asymptotic expansion of the solution are given in Section 5. Incidentally, 
an H-function representation of a one-sided Levy stable density is also obtained. 



2 Unified fractional kinetic equations 

In this Section, we present a method based on Laplace transform for deriving 
the solution of the unified fractional kinetic equations. 

Theorem 1. If Re(fj) > 0, dj > 0, j 6 N, and f(t) be a given function, defined 
on then the equation 

n 

N(t) - No f{t) = -J2 a i ° D t Ui N(t), (1) 
is solvable and its particular solution is given by 

- "ob-d 1 e (ri) , (0 ( ! r , {nW4 

1=0 n + ...+r„_i=i K r l)----K r n-1)- l M=1 J 

/"*/(«)(* [-ai(t-«r]d«, (2) 
Jo ^'Z^=i 

where the summation in (2) is taken over all nonnegative integers r\ , . . . , r„ 
such that n + . . . + r„_i = Z, and provided that the series and integral in (2) 
are converg cnt. Here D t 3 , j G N are Riemann-Liouville fractional integrals, 
defined by 

oDrf{t) = -f- /*(* - uy-\f(u)du, Re{v) > 0, (3) 
1 \ v ) Jo 

with D$f(t) = f(t) (Oldham and Spanier, 1974; Miller and Ross, 1993), 
Ep (z) is the generalized Mittag-Leffler function, defined by Prabhakar (1971) 
in terms of series representation as 

E f>J z ) = E T (3r + j)(r)[ (/3 ' 7 ' 6 G °' Re{(3) > °' Re{l) > 0) ' (4) 
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Proof. By the application of the convolution theorem of the Laplace transform 
(Erdelyi et al., 1953, p. 259) to (3), we find that 

L{ D^f(t);s} = 

= s-T{8), (5) 

where /~(s) = J °° e~ st f(t)dt,s e C,Re(s) > 0. Applying Laplace transform 
to (1) and using (5), it gives 

iV~( s ) = NoTM (6) 

1 + a-LS-"! + ... + a n s-^ 



\-^n—l —v- i i 



1 + ais~ Vl 



< 1. 



1=0 

If we employ the identity (Abramowitz and Stegun, 1968, p. 823) 

(*i + -+^)'= E (ri / )! (r y fW' w 

where the summation is taken over all nonnegative integers, ri,...,r n , such 
that n + . . . + r n = I, then for \ais~ Vl \ < 1, (7) transforms into the form 

it {s) = iVor( S )E(-i) f E (ri) , 1 l} , 

r 1 >...r n _ 1 >0 

(l + oia-^)'+ 1 ■ lj 

Taking the inverse Laplace transform of (8) by making use of the formula (Kil- 
bas, Saigo, and Saxena, 2004, eq. (12)) 

L- 1 {«-T(l - as- )- 6 ; t] = f- 1 ^^), (9) 

where i?e(s) > |a| 1 /«e(7) j fl e ( 7 ) > o,i?e(s) > 0, 

and applying the convolution theorem of the Laplace transform, the result (2) 
is established. 

Remark 1. The generalized Mittag-Lcffler function defined by (4) is studied by 
Prabhakar (1971) and Kilbas, Saigo, and Saxena (2004). Recently this function 
is used in the theory of finite-size scaling of systems with strong anisotropy and 
long-range interaction by Chamati and Tonshev (2006). 
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3 Special cases 

Some special cases of Theorem 1 are of interest to be highlighted. If we set 

v j = 3 v i a j = fflc^U £ N), we obtain 

Theorem 2. If Re(v) > 0, c > and f(x) G then the equation 

N( t ) N f( t ) = f2(rV rD t vr m), (io) 

r=l 

is solvable and its solution has the form 

N(t) = N j t J* f(u)E^[-c^(t - uY]du, (11) 

where E™ 1 (x) is the generalized Mittag-Lcfflcr function defined by (4) and pro- 
vided that the integral (11) is convergent. 

When n — 1, we obtain the following result given by Hille and Tamarkin 
(1930). 

Corollary 2.1. Let Re(v) > 0, c > and let f(x) G 9?+, then for the solution 
of the integral equation 

N(t) - N Q f(t) - -c v Q D^N(t), (12) 

holds the following formula 

N(t) = N j f J f{u)E v [-c v {t - u) v ]du, (13) 

where E u {z) is an entire function of order p = i and type a = 1, defined by 

E »W = E rfr^ + i) ' [v e C ' Re{v) > 0) ' (14) 
^i— i 

Note 1. The above result has also been given by the authors in a different form 
(Saxena, Mathai, and Haubold, 2004a, 2004b). 

If wc set f(t) = t'T- 1 Ei^[-(ct) v ], Theorem 2 yields 

Corollary 2.2. Let Re(v) > 0,Re(-f) > 0,c > 0, then for the solution of the 
kinetic equation 

n 

N(t) - N^El^-icty] = -J2(rV v oDt rv n G N (15) 
holds the relation 

N(t) = N^^El+^-tctYlne N. (16) 
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For f(t)=t p 1 , Theorem 2 yields the following result 

Corollary 2.3. If Re(p) > 0,Re(v) > 0,c > 0, then for the solution of the 
equation 

n 

N(t) - Not"- 1 = - Y,(rY V oDt r "N(t), r G N, (17) 

r=l 

holds the relation 

N(t) = N tP- 1 r(p)El\ p [-(ctnr G N. (18) 

For n = 1, cq. (18) reduces to a result given by Saxena, Mathai, and Haubold 
(2002, p. 283, eq. (15)). When a, = a j s uj , for j = 1, . . . , n, we obtain 

Theorem 3. Let Re{v) > 0, a > 0, t > 0, n > 1, |a™+ 1 s-( Il+1 ) 1/ < 1, and f(x) 
be a given function defined on 3? + , then the equation 

n 

N(t) - N f(t) = oD^N®, (19) 

r=l 

is solvable and its solution is given by 

jv(t) = iVo|^£/Hi? ( „ +1) ^K +1 (t- u )( n+1 ^]d u 

-a J (t - uy^E^^a^it - uf n+1 >}du} , (20) 
where E^ n+1 ^ u (z) is the generalized Mittag-Leffler function E a .p{z) defined as 

oc n 

E a , (z) = ]T Tfal + pY ( "' P G Cl Re{a) > °' Re{P) > 0) (21) 

and provided that the integral in (20) is convergent. 

If we take v 3 ' — jv, for j = 1, . . . n, then it is interesting to note that Theorem 
1 yields the following result given by (Miller and Ross, 1993) in a different form: 

Theorem 4. Let Re(v) > 0,aj > 0, and f(x) be a given function defined on 
5R+, |ais _l/ | < 1, then the fractional kinetic equation 

n 

N(t) - No f{t) = oD^Nit), (22) 

3=1 

is solvable and has the solution given by 

- £(-!)' E (ri) , (0 ( r lV (n( fl M + i) r "} (23) 
1=0 ri+ ... +rn _ 1= i ^ij.-..(r„_ij. j 

provided that the series and integral in (23) are convergent. 
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4 Fractional diffusion equation 

In this Section we present an alternative shorter method for deriving the solution 
of a diffusion equation discussed earlier by Kochubei (1990). 

Theorem 5. Consider the Cauchy problem 



with 



D?N(x, t) = -c"AN(x, t), (0 < a < 1; x € 5T; < t < T), (24) 



N(x,t = 0) = S(x),x e n, lim N(x,t) = (25) 

x — >±oo 



o-Df is the regularized Caputo (1969) partial fractional derivative with respect 
to t, defined by 



vD?N(x,t) 



1 



r(l-a) 



d f* N(x,s)ds N(x,0) 



dt 



(t - x) c 



t° 



and A is the Laplacian. The fundamental solution of the above Cauchy problem 
is given by 



N(x,t) = \x\- n Tt-*H{ 



2.0 



i(i,«) 



Ac u l(n/2,l),(l,l) 



(26) 



where H 1 ' 2 (.) is the H-function (Mathai and Saxena, 1978). 

Proof. Applying the Laplace transform with respect to t, using the result 
(Caputo, 1969) 

m— 1 

L{ D?N(x,t)} = s a N~(x,s) - s a " r " 1 ^ (r) (^-0), 

r=Q 

m — 1 < a < m, mGN, 
and Fourier transform with respect to x, gives 

s a N~* (k, s) - s"- 1 = -c v \k\ 2 N~' (jfe, s), 

where the symbol " ~ " indicates the Laplace transform with respect to the time 
variable t and the symbol "*" the Fourier transform with respect to the space 
variable x. 

Solving for N~ , we have 



„a-l 



iV~ (k,8) 



s a + c»\k\ 2 ' 



(27) 



By virtue of the following Fourier transform formula (Samko, Kilbas, and Marichev, 
1990, p. 538, eq. (27.1)) 
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(r. \ n/2 

(28) 

where the multidimensional Fourier transform with respect to x G 3?" is defined 
by 

{F x N)(t, t)= f N(x, t)e ixT dx (t G 5ft"; t > 0) (29) 
and is the modified Bessel function of the second kind, yields 



N(x, s) = C - V s a - X 



(27T)" 



„a/2 



n-2/2 



S 2 x 



(30) 



In order to invert the Laplace transform, we employ the following result 
given by the authors (Saxena, Mathai, and Haubold, 2006) 



L- 1 { 8 -"K v (z8"y,t} = -t p - l Hl% 



zH~ 2 ° 



(P,2<r) 

i) 



(31) 



where K v (x) is the modified Bessel function of the second kind, Re(z 2 ) > 
0, Re(s) > 0. Thus we obtain the solution in the form 



N(x,t) 


= 2^ 


) 2 C 2 


4 X 1 2 t 2 


y. n 
4 




tt2,0 
H l,2 


"i-«|a;| 2 


(l-f-^,a) 


(32) 




4c" 


(^,1),(^,1). 


By virtue of a result in Mathai and Saxena (1978), 






|(o p ,Op) 


rrm,n 


:(a p +aA p ,A p ) 
X \(b q +aB q ,B q ) _ 


(33) 


the power of the expression [ji " 


a;| 2 } /4c" 


] can be absorbed inside the H- 



function and consequently we obtain 

N(x,t) = \^x\- n Hl>° 2 
Remark 1. If we employ the identity (Mathai and Saxena, 1978) 



~t- a \x\ 2 


(La) 


4c" 


(f,l),(l,l). 



(34) 



Tjm.n 



\,(a p ,A p ) 



Tjm.n 



,(a p ,A p /X) 
kb q ,B q /\) 



,A > 



the solution given by (32) can be expressed in the form 



1 



N{x,t) = -iTT'd - " Ht% 
a 



2 \ 1 / a 



(i,i) 



(35) 



(36) 
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where a > 0. We also note that the above form of the solution is due to 
Schneider and Wyss (1989). There is one importance of our result (32) that it 
includes the Levy stable density in terms of the H- function as shown in (45). 
Similarly, using the identity (35) we arrive at 

N{c,t) = \\^x\- n Hl>° 2 

where n is not an even integer. This form of the H-function is useful in deter- 
mining its expansion in powers of x. Due to importance of the solution, we also 
discuss its series representation and behavior. 



t-~\x\ 


(i,f) 


2c% 





(37) 



5 Series representation of the solution 

Using the series expansion for the H-function given in Mathai and Saxena (1978), 
it follows that 



rr2,0 
H l,2 



(1,1) 

(9.£).(i.£) 



2m J L T(l-s) XdS (38) 



^ r(i-f -0(-i)V*(t+o ~ r(% - i -Q(-i)'x"( 1+ ') 

h r(l-^ -«/)(/)! + h r(l-a-a/)(Z!) 



where n is not an even integer. 
Thus for n = 1, we find that 

1 °° A" 

where A — ^ and the duplication formula for the gamma function is used. 

For n = 2, the H-function of (37) is singular and in this case, the result is 
explicitly given by Barkai (2001) in the form 

N(x,t)~ rn 1 ,.> [—]. (40) 
7rl (1 — a)r" x 

For n = 3, the series expansion is given by 

N{xt) = I f (-V lAl/2 (41) 

{ , ) 4tt^/ 2 A1/ 2 Z-r[l-a(l + //2)]" 1 J 

From above, it readily follows that for n — 3 and ct^ 1 

JV(a:,i) ~-, as noo. (42) 

a; 
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It will not be out of place to mention that the one sided Levy stable density 
can be obtained from Laplace inversion formula (31) by virtue of the identity 



(43) 



and can be conveniently expressed in terms of the Laplace transform 



Jo 



/ e- ut § p {t)dt = e~ 



Re(u) > 0, Re(p) > 0. 



(44) 



The result is 



M*) = -His 



(p > 0). 



(45) 



Note 2. This result is obtained earlier by Schneider and Wyss (1989) by fol- 
lowing a different procedure. Asymptotic behavior of $ a (t) is also given by 
Schneider (1986). 

In conclusion, we mention that some of the results derived in this article 
may find some applications in problems associated with models of long-memory 
processes driven by Levy noise and other related problems, see the article by 
Anh, Heyde, and Leonenko (2002). 
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